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Abstract 

Accumulation point of period-tripling bifurcations for complexified Henon map is found. 
Universal scaling properties of parameter space and Fourier spectrum intrinsic to this critical 
point is demonstrated. 

It is known that complexification of real ID logistic map 

z n+l — A — Z„ 



■'n i 



where A,z6C leads to the origination of the Mandelbrot set at the complex parameter A plane pQ and 
a number of other accompanying phenomena. So, except for the transition to chaos through cascade 
of period-doubling bifurcations, intrinsic to the real logistic map, other bifurcation sequences are 
possible. For example, the universal properties of transition to chaos through period-tripling cascade 
is studied in paper of Golberg-Sinai-Khanin j2J. Corresponded critical point (GSK point) is situated 
at 

A c = 0.0236411685377 + 0.7836606508052i (2) 

and characterized by following critical indexes namely, by critical multiplier // c , scale factor a 
and parameter scaling constant 5 (see Fig. 1(a)): 

p c = -0.47653179 - 1.05480867?, (3) 

a = -2.0969 + 2.3583i, (4) 
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.5 = 4.6002 - 8.9812i. (5) 

Opportunity for realization of the phenomena, characteristic for the dynamics of complex maps 
(Mandelbrot set etc.) at the physical systems seems to be interesting problem jH]. In the context of 
this problem the following question is important: Does phenomena of dynamics of the ID complex 
maps like classic Mandelbrot map survive for the two-dimensional maps. For example, from the 
point of view of possible physical applications, more realistic model rather than logistic map, is the 
Henon map 

z n +i = f{z n , w n ) = A - z 2 n -d- w n , w n+l = g(z n , w n ) = z n . (6) 

In the real variable case the system © is 2D invertible map and, hence, can be realized as Poincare 
cross-section of flow system with three dimensional phase space - minimal dimension, providing 
opportunity of nontrivial dynamics and chaos. Henon map is suitable for modelling of the chaotic 
dynamics of the generator with non-inertial nonlinearity, dissipative oscillator and rotator with peri- 
odic impulse driving force etc. |4 . Moreover, Henon map expresses the principal properties of large 
class of differential systems. 

Let us complexify the map (0) in a such way that z, w, A G C, d e R. According to the work j3Ej, 
complexified Henon map can be reduced to the two symmetrically coupled real Henon maps and can 
be realised at the physical experiment [Zj. 

Let us remark that with \d\ < 1 the Henon map is dissipative system, with |d| — > 1 - it is area- 
preserving map, and with d — > it corresponds to complex ID quadratic map (0), which describes 
the universal scenario of transition to chaos through period-tripling bifurcations. In present work we 
aim to be convinced of existence of GSK critical point for the Henon map with 1 > \d\ ^ 0, to find 
this point and to demonstrate intrinsic to it scaling characteristics. 

Let us shortly explain essence of the numerical procedure for the calculation of the GSK point. 
From the work jH] , the simple and effective method of calculation of the critical point of transition 
to chaos though period-doubling bifurcations for ID maps and also for the more general systems is 
known. We have generalized this method to the case of GSK critical point and have applied it to the 
2D complex Henon map. The method is based on the fact, that at the GSK point the infinite number 
of the unstable cycles with tripling periods must exist, and multipliers of these cycles must tend to 
the universal constant /i c . Therefore, for the finding of the point Xf^° it is enough to calculate with 
fixed parameter d the point at the A plane, satisfying to these conditions. Therefor it is necessary to 
solve the system of nonlinear complex equations: 

f N (z N ,w N ) = z N , g N (z N ,w N ) = w N , ii(z n , w N ) = fM c (7) 
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concerning the cycle elements zjq, w^, where N = 3 k (k — > oo) and parameter A. The system of 
equations (JZJ) can be solved numerically by Newton procedure. The multiplier /j,(zn,wn) can be 
found as maximal by absolute value eigenvalue of the Jacobi matrix J(zn,wn), where 

( df N (z,w) dg N (z,w) \ 
of (z,w) ag (z,w) j 
dw dw / 

Let us mark, that for the best convergence of a Newton method, expediently to use inheriting 
starting conditions. At the each consequent step of calculations, as a starting position of the critical 
point the value of the parameter A obtained at the previous step is used. As starting position of 
the elements of cycles the values obtained according to the universal scaling properties of the phase 
space at the critical GSK point are used. First elements of different cycles selected properly should 
obey the equations 

zn = z c + C\a~ N , wn = w c + C^a ', (9) 

where z c and w c - are the constants, which correspond to the scaling center at the phase space, and 
C±, C*2 -are some complex coefficients. In the case of the ID complex map, i.e. with d = 0, the 
scaling center is situated at the origin. With d 7^ the scaling center appears displaced concerning 
an origin of coordinates and can be found from the following conditions 

z c = (az N+1 - z N )/(a - 1), w c = (aw N+1 - w N )/(a - 1). (10) 

Let us summarize the algorithm of numerical calculations. At the first stage, one must define 
the critical GSK point for the cycles with not very large periods and dragging the parameter d from 
(critical point of the ID map) to the necessary value. This allow to use the known critical GSK 
point (J2J) as a starting condition at the beginning of the calculation procedure at d = 0. Then, with 
fixed parameter d one must consequently improve the critical point using the more long-periodic 
cycles. 

Numerical calculations have shown that GSK point for the Henon map © with traditional, 
mentioned in the original work of Henon 9j value of parameter d = —0.3 is situated at 

A d=-o.s = _o.24388583757 + 0.69478896727i. (11) 

The position of the scaling center is defined by the following values 

Zc = -0.05289397 + 0.15014466i, w c = 0.56676843 + 0.72274977i. (12) 

In the neighborhood of the point (fTTj) at the parameter plane (ReA, ImA) the scaling, characterized 
by constant 5 takes place (see Fig. 1(b)). The existence of the infinite number of the cycles with 



tripling periods immediately at the critical point Af =_0 ' 3 is proved to be true by the structure of 
the Fourier spectrum. At Figure 2 the dependence of the spectral intensities versus the logarithm 
of frequency / is represented, and amplitudes of harmonics are rated to the f K , where k — 6.15 - 
is universal constant, characterizing the self-similar fractal structure of the spectrum at the GSK 
point [TO] . 
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Figure 1: Demonstration of the scaling properties of the Mandelbrot set near the GSK point (critical 
point is situated at the center of the small fragments) for the complexified Henon map (5) at the 
plane of complex parameter A with d = —0.0 (a) and d = —0.3 (b). Performed by the multiplication 
to the complex scaling constant 5 the fragment of the Mandelbrot set has the same structure with 
the previous fragment. 
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Figure 2: Double logarithmic plot of Fourier spectrum of the signal, originated by the map (jBJ) at 
the critical point pip. Harmonics of the tripling frequencies has the same amplitudes. 
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